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摘要 : 2020 年 7 月 31 日 上 午 ， 北 斗 三 号 全 球 卫 星 导 航 系统 正式 开通 。 随 着 北斗 系统 的 逐步 完 
善 以 及 精度 的 不 断 提 升 ， 用 户 对 北斗 系统 的 完好 性 要 求 也 越 来 越 高 。 北 斗 三 号 系统 在 广播 星 历 


中 利用 SISAls, SISA c, SISAIo 和 95719471oc 这 4 个 参数 表示 星 历 和 钟 差 参 数 精度 ， 但 
是 目前 4 个 SISAI 参数 尚未 在 接口 控制 文件 (Interface Control Document, ICD) 中 给 出 用 户 


H 


算法 ， 使 北斗 三 号 系统 的 完好 性 服务 不 能 应 用 于 实际 。 研 究 了 北斗 三 号 系统 基本 完好 性 参数 
SISAI 的 等 级 转换 算法 和 用 户 使 用 算法 ， 并 根据 国际 民航 组 织 (International Civil Aviation 
C Organization, ICAO) 非 精 密 进 近 (Non-Precision Approach, ICNPA) 阶段 完好 性 要 求 ， 利 用 国 
际 GNSS 监测 评估 系统 (International GNSS Monitoring and Assessment System, IGMAS) 全 
球 网 实测 数据 和 北斗 三 号 卫星 广播 星 历 中 SISA 参数 ， 对 北斗 三 号 系统 的 完好 性 服务 性 能 进行 
了 初步 验证 分 析 。 试 验 结果 表明 ， 随 着 北斗 三 号 系统 组 网 完成 ， 利 用 SISA 参数 转换 和 用 户 使 
CN 用 算法 ， 北 斗 三 号 系统 NPA 阶段 完好 性 服务 能 力 为 100%。 
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全 球 卫 星 导 航 系 统 (Global Navigation Satellite System, GNSS) 服务 性 能 的 主要 指标 包 
括 精度 、 完 好 性 、 连 续 性 和 可 用 性 。 其 中 导航 系统 完好 性 是 指 在 导航 系统 提供 的 导航 信息 不 
能 满足 服务 性 能 要 求 时 (导航 系统 发 生 故 障 或 者 误差 超 限 ， 无 法 用 于 导航 定位 )， 给 用 户 提 
供 及 时 告警 的 能 力 。 完 好 性 是 衡量 卫星 导航 系统 是 否 可 靠 的 一 个 重要 指标 ， 随 着 全 球 卫 星 导 
航 系统 服务 精度 的 不 断 提升 ， 应 用 范围 不 断 扩 大 ， 关 系 用 户 使 用 安全 的 系统 完好 性 问题 ， 越 
来 越 受 到 重视 ， 各 GNSS 系统 都 将 系统 完好 性 建设 作为 系统 升级 的 重要 工作 之 一 “。 
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卫星 导航 系统 的 完好 性 监测 手段 主要 包括 卫星 自主 完好 性 监测 (Satellite Autonomous 
Integrity Monitoring, SAIM)、 地 面 完 好 性 监测 (GNSS Integrity Channel, GIC) 和 接收 机 自 
主 完好 性 监测 (Receiver Autonomous Integrity Monitoring, RAIM)。 其 中 ， 卫 星 导 航 系统 
的 建设 方 主要 依靠 GIC 和 SAIM 技术 实现 系统 的 完好 性 监测 ， 而 RAIM 技术 则 通过 用 户 接 
收 机 算法 完成 对 故障 卫星 的 识别 。GIC 技术 是 指 在 利用 系统 布设 的 地 面 监测 站 ， 收 集 导 航 
卫星 的 观测 数据 ， 对 导航 电文 的 空间 信号 精度 进行 预测 和 监测 ， 生 成 相应 的 完好 性 电文 ， 并 
随 导 航 电 文 一 起 播发 给 用 户 ， 用 户 利 用 完好 性 电文 参数 对 卫星 健康 状态 和 导航 系统 播发 信 
息 的 可 靠 性 进行 判断 ”“。SAIM 技术 则 是 在 导航 卫星 上 直接 对 本 星 信号 的 异常 进行 监测 ， 
根据 监测 结果 ， 播 发 卫星 信号 健康 标识 ， 电 文 参数 简单 ， 告 警 实 时 性 高 ， 是 下 一 代 卫 星 导 航 
系统 完好 性 监测 的 重要 手段 ”。RAIM 技术 是 接收 机 在 利用 所 接收 的 多 颗 导 航 卫 星 观 测量 
进行 定位 处 理 的 同时 ， 采 用 最 小 二 乘法 或 奇偶 空间 矢量 算法 原理 对 多 星 元 余 观 测量 进行 故 
障 卫星 的 检测 和 排除 分 析 处 理 ”， 该 技术 一 般 只 能 针对 单 卫星 较 大 故障 进行 监测 ， 而 且 会 降 
低 系 统 可 用 性 ， 因 此 完好 性 监测 能 力 相 对 较 弱 。 

BDS, GPS, Galileo 等 GNSS 导航 系统 向 用 户 播发 的 完好 性 电文 参数 中 ， 一 类 重要 信息 
为 广播 星 历 的 卫星 轨道 和 钟 差 的 预测 精度 信息 。 由 于 监测 体制 和 电文 设计 不 同 ， 各 GNSS 
系统 播发 的 完好 性 电文 参数 并 不 相同 。 早 期 系统 实现 ， 通 常用 一 个 完好 性 电文 参数 表示 卫 
轨道 和 卫星 钟 差 的 综合 径 向 预测 精度 ， 随 着 系统 完好 性 监测 能 力 提 升 ， 进 一 步 考 虑 卫星 轨 
和 卫星 钟 差 的 不 同 误差 特性 ， 采 用 不 同 的 完好 性 参数 分 别 反映 卫星 轨道 和 卫星 钟 差 的 电 
精度 。 其 中 GPS 系统 播发 导航 电文 中 利用 URAT (User Range Accuracy Index) 数 反映 广播 
星 历 空间 信号 预测 精度 ， 在 新 设计 播发 的 CNAV (Civil Navigation) 电文 中 ， 利 用 URATa， 
DVR47neao U RAlscai, U RAInca2 这 4 个 参数 ， 分 别 表 示 星 历 切 平 面 预报 误差 、 星 历 径 向 与 
钟 偏 预报 误差 、 钟 速 预报 误差 和 钟 漂 预 报 误差 ” ;Galileo 系统 和 BDS-2 (北斗 区 域 系统 ) 系 
统 分 别 利 用 SISAI (Signal In Space Accuracy Index) 参数 ”和 URAI 表示 广播 星 历 空间 
信号 预测 精度 ”" ; BDS-3 (北斗 全 球 系统 ) 系统 对 完好 性 监测 体制 进行 了 全 面 升 级 ， 利 用 
SISAI,,, SISAL, SISAL, SISA o 这 4 个 参数 表示 星 历 切 平面 预报 误差 、 星 历 径 向 
与 钟 偏 预 报 误差 、 钟 速 预报 误差 和 钟 漂 预 报 误 差 ” 7. $ GNSS 系统 在 广播 电文 中 播发 的 
星 历 和 钟 差 完好 性 电文 参数 和 含义 如 表 1 所 示 。 


KOodk mu 


© R1 GNSS 系统 星 历 与 钟 差 完 好 性 电文 参数 对 比 
系统 完好 性 参数 含义 
GPS (LNAV) URAI 用 户 距 离 精度 指数 
URAR 道 平面 预测 精度 指 妆 
GPS (CNAV) d REDEUNT NES 
URAlIca0, U RAlsca1, U RAIncaa 卫星 钟 差 预测 精度 指数 
Galileo SISAI 空间 信和 号 精度 指数 
BDS-2 URAI 用 户 距 离 精度 指数 
—€— SISAI,。 卫星 轨道 的 切 向 与 法 向 精度 指数 
E 3S7T547cb， 卫 星 轨 道 的 径 向 及 卫星 钟 固 定 偏差 
97947ocb， 97947ocl， SISAloa 精度 指数 ; SISAlsa:; 1 星 钟 频 ji 精度 ES 
S7T5S47ca， 卫 星 钟 频 漂 精度 指数 


m= 
p= — 
— 
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系统 完好 性 的 服务 性 能 一 般 以 告警 限 值 、 告 警 时 间 和 和 危险 误导 信息 (Hazardously 
Misleading Information, HMI) 概率 具体 表示 。 告 警 限 值 是 指 对 应 飞行 阶段 保证 安全 操作 的 
定位 误差 限 值 ， 分 为 水 平 告警 限 值 (Horizontal Alert Limit, HAL) 和 垂直 告警 限 值 (Vertical 
Alert Limit, VAL) 两 类 。 告 警 时 间 是 指 从 系统 出 现 故 障 开始 到 用 户 收 到 告警 所 允许 的 最 大 
WEJ. HMI 概率 是 指 当前 定位 值 落 在 告警 限 值 范围 之 外 的 飞行 危险 概率 ， 一 般 要 求 不 
大 于 10-7h-1”” 。 由 于 涉及 生命 安全 服务 ， 航 空 用 户 对 卫星 导航 系统 完好 性 提出 了 最 严 
格 的 要 求 ， 表 2 给 出 国际 民航 组 织 对 民航 各 飞行 阶段 的 导航 完好 性 性 能 要 求 。 


表 2 ”国际 民航 各 阶段 导航 完好 性 性 能 需求 表 


m 完好 性 
飞行 阶段 — T — 
完好 性 风险 告警 限 值 /m 告警 时 间 /s 
非 精 密 进 近 (NPA) 107757! 556 (H AL) 10 
有 垂直 引导 进 近 (APV-I) | 2x 10-" 每 次 进 近 556 (HAL); 50 (VAL) 10 
有 垂直 引导 进 近 (APV- ID) | 2x 10- 每 次 进 近 40 (HAL); 20 (VAL) 6 
CAT I 精 密 进 近 2 x 10-7 每 次 进 近 40 (HAL); 20 (VAL) 6 


如 表 1 所 示 ， 各 GNSS 系统 通过 广播 星 历 播发 的 完好 性 电文 参数 均 是 索引 值 ， 用 户 需 
要 根据 GNSS 系统 发 布 的 用 户 接口 控制 文件 将 完好 性 参数 索引 值 转换 为 具体 的 空间 信和 号 精 
度 预 测 值 ， 并 根据 用 户 算法 进行 应 用 。BDS-3 系统 在 2018 年 12 月 27 日 开通 了 基于 18 颗 
MEO 卫星 的 基本 服务 ， 于 2020 年 7 月 31 日 正式 开通 基于 24MEO/3IGSO/3GEO 星座 的 全 
球 服务 。 目 前 BDS-3 系统 尚未 在 其 ICD 中 给 出 各 SISAI 索引 参数 的 等 级 转换 关系 ”， 也 
未 给 出 SISA 参数 的 用 户 应 用 算法 。 本 论文 设计 了 BDS-3 系统 SISAI 参数 的 等 级 转换 算 
法 及 用 户 应 用 算法 ， 基 于 定位 域 ， 研 究 了 SISA 参数 的 完好 性 评估 方法 。 利 用 IGMAS 全 
球 监测 网 实测 数据 ， 对 比 不 同 阶段 BDS-3 系统 完好 性 服务 能 力 ， 首 次 对 BDS-3 系统 播发 的 
SISA 参数 的 完好 性 性 能 进行 了 初步 验证 。 


一 


2 北斗 三 号 系统 SISAL 参数 索引 转换 算法 


受到 导航 电文 播发 带宽 的 限制 ，GNSS 系统 发 给 用 户 的 空间 信号 完好 性 参数 均 以 等 级 
索引 值 的 形式 表示 。 在 BDS-3 系统 广播 星 历 中 ， 利 用 完好 性 参数 SISAL 表达 广播 星 历 中 
卫星 轨道 和 钟 差 的 预报 精度 。 卫 星 通过 广播 电文 播发 的 完好 性 参数 索引 值 包括 : SISAL., 
SISAl«4, 515471uca。 用 户 在 使 用 完好 性 参数 前 ， 需 要 将 完好 性 参数 的 索引 值 转化 为 具体 的 
预测 精度 值 ， 目 前 北斗 ICD 中 尚未 发 布 各 个 参数 专门 的 等 级 对 应 关系 表 。 
因 卫 星 轨 道 与 钟 差 误差 的 特性 不 同 ， 对 不 同位 置 用 户 的 影响 不 同 ， 其 中 卫星 轨道 切 向 
与 法 向 方向 预报 精度 以 SISA 参数 表达 ， 卫 星 轨 道 径 向 及 钟 差 预报 精度 以 SISA 
参数 表达 。 同 时 ， 根 据 钟 差 误差 的 二 次 模型 特性 ， 又 将 SISAI 参数 进一步 细 化 为 
SISA, SISAl«4 和 SISA. SISAL 参数 综合 表示 卫星 钟 偏差 和 卫星 轨道 径 癌 
RÆ, SISAhba 表示 卫星 钟 钟 速 误差 ，STSATcs 参数 表示 卫星 钟 漂移 误差 。 
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其 中 SISA. R 
SISA» 索引 值 范 围 
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一 16 ~ 十 15， 不 同 索 引 值 表示 的 参数 范围 如 表 3 所 示 。 


也 为 -16 ~ 十 15， 不 同 索 引 值 表 示 的 参数 范围 SISA. 参数 完全 相 


H, ABER. SISAlL,A4 和 SISAL 参数 范围 为 0 ~ 7， 表 示 8 个 等 级 ， 没 有 专门 的 等 级 
对 应 关系 表 。 
表 3 SISA. 索引 关系 表 
SISAI 索引 SISAoe/m SISAI。。 索引 SISAoe/m 

15 6 144.00 < SIS Aoe 一 1 1.20 < SISAoe < 1.70 
14 3 072.00 < SIS As. «6 144.00 三 之 0.85 < SIS Ave < 1.20 
13 1536.00 < SISA,, < 3 072.00 一 3 0.60 < SISA,. < 0.85 
12 768.00 < SIS Ave «1 536.00 一 4 0.43 < SIS Ave < 0.60 
11 384.00 < SIS Ave «768.00 一 5 0.30 < SIS Ave < 0.43 
10 192.00 < SIS Aoe «384.00 —6 0.21 < SIS Ave < 0.30 
9 96.00 < SIS Ave «192.00 一 7 0.15 < SIS Ave < 0.21 
8 48.00 < SIS As. < 96.00 一 8 0.11 < SIS Aoe < 0.15 
7 24.00 < SIS Ave < 48.00 一 9 0.08 < SIS Aoe < 0.11 
6 13.65 < SIS Aoe < 24.00 —10 0.06 < SIS Ave < 0.08 
5 9.65 < SIS Ae. < 13.65 —11 0.04 < SIS Ave < 0.06 
4 6.85 < SIS Aoe < 9.65 一 12 0.03 < SIS Ave < 0.04 
3 4.85« SIS Ae. < 6.85 —13 0.02 < SIS Ave < 0.03 
2 3.40 < SIS Aoe < 4.85 一 14 0.01 < SIS Aoe < 0.02 
1 2.40 < SIS Ae. < 3.40 一 15 0< SISAoe < 0.01 

0 1.70 < SIS Aoe «2.40 一 16 没有 精度 预测 ， 使 用 有 风险 


对 于 SISAL. $8 SISAl; 参数 ， 本 文 设计 按照 索引 值 等 级 范围 内 取 标 准 优 
体 转换 如 式 (1) 所 示 : 
SIS Aoejocb = 2(11N/2) ， 

SIS Aoejocb = 209, 

SIS Aoc/ocb = 2.8, 5.7, 11.3 ， 


其 中 ， 为 表 3 中 的 索引 参数 值 ， 式 (1) 中 N 小 于 0 时 取 对 应 表 3 中 各 索引 值 对 应 精度 的 


引 值 对 应 的 等 级 转换 值 


ERRE. 


(N < 6) 
(N>6, N<15) , 
(N =1, 3, 5) 


对 于 SISAI4 M SISA 参数 ， 索 引 值 的 等 级 值 转换 如 式 (2) 所 示 : 


利用 式 (1) 和 式 (2)， 可 将 导航 电文 


示 值 。 


SIS Aoc, = 2- (8184113 
SIS AS = 27 (8184121 


PERIE SISAI 等 级 索引 值 转换 为 


3 SISA 参数 用 户 应 用 算法 


用 户 在 收 到 广播 星 历 中 完好 性 


FE 参数 的 索引 值 后 ， 首 


H14) 


-28) 


先 将 索引 值 转换 为 具体 的 精度 值 。 
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对 于 卫星 钟 差 的 预报 误差 ， 由 3 个 参数 表示 ， 需 要 按照 钟 差 预报 模型 ， 利 用 索引 转换 后 的 
SIS Aoc, 5754。 和 SIS Aoo 参数 计算 卫星 钟 差 的 预报 精度 ， 如 式 (3) 所 示 : 


| SIS Aoc = SIS Aoc + SISA«a(t— top), t—top < 93600 


SIS Aoc = SIS Aoc + SIS Aoc lt — top)+ ; (3) 
SIS Aoc2lt — top —93600)?, — t— top > 93600 


其 中 ，t 为 系统 时 间 ， 刀 p 为 电文 参考 时 间 ， 单 位 都 是 秒 。 

用 户 获得 95154。。 和 SISA 参数 后 ， 还 需要 根据 误差 特性 计算 SISA 综合 参数 ， 以 反 
映 广 播 星 历 的 空间 信号 精度 。 根 据 轨道 平面 误差 和 轨道 径 向 误差 、 钟 差 误差 在 用 户 最 差 位 置 
方向 的 近似 投影 ， 采 用 式 (4) 表示 SISA 综合 参数 : 


SISA = J(SISA,, x sin14)?+ SISA2, . (4) 


用 户 利用 系统 播发 的 完好 性 参数 ， 计 算 定位 保护 限 值 ， 以 反映 根据 观测 量 解 算出 来 的 位 
置 结果 的 可 靠 性 。 系 统 的 完好 性 服务 性 能 可 通过 定位 误差 保护 限 值 (Protection Level, PL). 
告警 限 值 (Alert Limit, AL) 和 定位 误差 (Position Error, PE) 这 3 个 参数 表达 。 定 位 误差 是 
指 用 户 利用 观测 数据 和 广播 星 历 计算 的 位 置 与 真实 位 置 的 偏差 ， 定 位 保护 限 值 是 指 用 户 利 
用 系统 播发 的 完好 性 参数 计算 的 定位 误差 的 保守 估计 ; 告警 限 值 是 指定 位 误差 的 告警 门限 ， 
若 计算 的 定位 保护 限 值 超 过 告警 限 值 ， 则 表示 导航 系统 的 服务 “不 可 用 ”。 告警 限 值 是 由 导 
航 定 位 服务 的 行业 来 确定 ， 目 前 国际 民航 根据 飞机 航路 的 不 同 飞行 阶段 ， 给 出 了 不 同 告警 限 
值 要 求 ， 如 表 2 所 示 。PE, PL 和 AL 的 关系 如 图 1 所 示 。 


保护 水 平 
(HPL/ VPL) 


真实 定位 误差 
(HPE/ VPE) 


告警 限 值 
(HAL/ VAL) 


图 1 PL/PE/AL 三 者 关系 


在 实际 情况 下 ， 用 户 (尤其 是 实时 动态 用 户 ) 通常 并 不 知道 自己 的 真实 位 置 ， 也 就 无 法 
计算 自己 的 导航 定位 误差 ， 只 能 通过 完好 性 参数 计算 相应 的 定位 保护 限 值 作为 定位 误差 的 
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保守 估计 。 而 卫星 导航 系统 需要 确保 ， 利 用 系统 播发 的 完好 性 参数 计算 的 保护 限 值 应 该 大 于 
j 户 真实 的 定位 误差 ， 若 利用 完好 性 参数 计算 的 保护 限 值 小 于 用 户 告 警 限 值 ， 但 用 户 真实 定 
立 误差 大 于 用 户 的 告警 限 值 ， 则 表示 发 生 了 完好 性 风险 事件 。 卫 星 导航 系统 最 重要 的 完好 性 
性 能 之 一 就 是 完好 性 风险 概率 。 

在 具体 应 用 中 ， 导 航 用 户 的 定位 误差 又 可 分 为 水 平定 位 误差 (Horizontal Position Error, 
HPE) 和 垂直 定位 误差 (Vertical Position Error, VPE)， 不 同 用 户 对 水 平和 垂直 方向 的 定位 
精度 要 求 不 同 ， 因 此 将 告警 限 值 分 为 水 平 告警 限 值 (HAL) 和 垂直 告警 限 值 (VAL)。 相 应 地 
用 户 利用 完好 性 参数 计算 的 保护 限 值 也 可 分 为 水 平 保护 限 值 (Horizontal Protection Level, 
HPL) 和 垂直 保护 限 值 (VPL)。 五 PL 和 VPL 的 具体 计算 公式 如 下 : 


ea 


d2 2 2 
HPL= EA east — | (e me 十 di ; (5) 
VPL= Ky: dy, (6) 


RP, Ky 和 Ko 分 别 为 根据 系统 完好 性 服务 等 级 的 水 平 风险 概率 Py umi 和 垂直 风险 概率 
Passi 得 到 的 置信 概率 相对 应 的 系数 ， 根 据 国 际 民航 公约 附件 10”， 飞 机 不 同 飞 行 阶段 的 
完好 性 风险 概率 系数 值 为 : Kunpa = 6.18, Ky pa = 6.0, Ky pa = 5.33. XX (5), (6) 中 di; 为 
利用 空间 几何 投影 矩阵 将 伤 距 域 视线 方向 残 差 投影 到 定位 域 的 结果 ”' ， 可 表示 为 


d?a den der id 


2 
dEN duh dyu dNT 


-(GTWGqG)'aG"W , (7) 
dgu duu dr 2 


dgr dyr dur dT 
式 中 ，G 为 观测 和 矩阵， 第 i 行 的 元 素 为 : 
G = [-cos(Ele)sin(Azi) —cos(Ele)cos(Azi) —sim(Ele) 1] . (8) 
BÓBEE W 是 对 角 阵 ， 可 表示 为 : 


1 1 
W - di " 9 
Tb (zc + SISA zn ssa] (9) 


Kb, SISA: 表示 各 卫星 的 S154 综合 参数 ， 由 式 (4) 计算 得 到 ，o, 为 传播 段 观测 误差 ” ， 
02,7 m 02i 3E T2 opo,i? 

其 中 ， 用 户 设 备 误差 car; 和 对 流 层 残余 误差 otropo,i 由 用 户 设备 和 对 流 层 模型 决定 ， 表 
示 如 下 : 


1 


x0.12 , (10) 


2 ze 
Ütropo,i m 


0.002 + sin?(0;) 
= 0.15 + rand(0.07) , (11) 


Oair, i 


其 中 ，4b 为 高 度 角 。 


-H 
H 


H 


F! 
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4 SISA 参数 完好 性 评估 方法 


BDS-3 系统 完好 性 服务 目标 为 满足 民用 航空 非 精密 进 近 NPA 阶段 的 完好 性 要 求 ，NPA 
服务 主要 关注 水 平方 向 的 完好 性 服务 能 力 。 利 用 IGMAS 全 球 分 布 监测 的 实测 数据 ， 按 照 斯 
坦 福 图 ”方法 对 北斗 系统 SISA 参数 的 完好 性 性 能 进行 验证 。 如 图 2 所 示 ， 按 照 斯 坦 福 回 
评估 方法 ， 将 系统 的 完好 性 服务 性 能 划分 为 四 种 情况 ”: 


风险 : AL<PL<PE 


正常 导航 : PE«PL«AL 


PE 


20] 


2 ”完好 性 斯 坦 福 分 析 图 ' 


(1) 当 定 位 误差 小 于 定位 保护 限 值 且 定位 保护 限 值 小 于 告警 限 值 时 ， 即 HPE < HPL < 
万 4 时， 表示 完好 性 服务 “可 用 ”， 由 图 2 中 白色 区 域 表示 。 

(2) 当 定位 误差 大 于 定位 保护 限 值 ， 但 定位 误差 小 于 告警 限 值 时 ， 即 HPL < HPE < 
及 AL 时， 表示 完好 性 服务 出 现 了 “ 漏 警 >， 由 图 2 中 蓝 色 区 域 表示 。 

(3) 当 定 位 误差 大 于 告警 限 值 ， 且 定位 保护 限 值 小 于 告警 限 值 时 ， 即 HPL < HAL < 
万 PE 时 ， 表 示 完 好 性 服务 出 现 了 “严重 的 完好 性 风险 ”， 图 2 中 红色 区 域 表 示 。 

(4) 当 定 位 保护 限 值 大 于 告警 限 值 时 ， 即 HPL > 及 4L， 表 示 导 航 服务 “不 可 用 ”， 由 
2 中 黄色 和 灰色 区 域 表示 。 其 中 黄色 区 域 表示 没有 完好 性 漏 警 ， 灰 色 区 域 表示 存在 完好 性 
漏 警 。 

按照 斯 坦 福 图 评估 方法 ， 选 择 测 站 坐标 已 知 的 观测 数据 进行 系统 完好 性 服务 的 评估 。 
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IGMAS 测 站 的 坐标 可 利用 精密 定位 等 方法 进行 精确 标定 ， 利 用 BDS-3 系统 播发 导航 电文 计 
算 各 测 站 的 位 置 ， 与 精确 坐标 进行 比较 ， 可 以 获得 各 测 站 的 导航 定位 误差 


POSerr = POSpost — POSnay , (12) 


式 中 ，POSerr 表示 定位 误差 ，POSbost 为 事后 精密 标定 的 测 站 坐标 ，POSuaav 表示 利用 导 
航 电文 计算 的 测 站 坐标 。 
利用 导航 电文 播发 的 S154 参数 ， 按 照 第 3 章 方法 计算 水 平 保护 限 值 五 PL， 分 别 对 完 
好 性 服务 的 可 用 性 、 漏 警 率 、 严 重 的 完好 性 风险 和 服务 不 可 用 的 概率 进行 统计 分 析 。 
BDS-3 系统 的 NPA 完好 性 服务 可 用 性 评估 方法 为 : 
ypa NHAL> HPL>I(POSerral 
Nan 
BDS-3 系统 的 NPA 完好 性 服务 漏 警 概率 分 析 ， 即 分 析 水 平定 位 误差 大 于 水 平 保护 限 
值 ， 但 小 于 水 平 告警 限 值 的 概率 : 
mr, — NHAL > |(POSerra)| > HPI) 
Nan 
BDS-3 系统 的 NPA 完好 性 服务 完好 性 风险 概率 分 析 ， 即 分 析 定 位 误差 大 于 水 平 告警 限 
值 ， 定 位 保护 限 值 小 于 水 平 告警 限 值 的 概率 : 
N(I(POSerraj| > HAL > HPL) 
Nan l 
BDS-3 系统 的 NPA 完好 性 服务 不 可 用 概率 分 析 ， 即 分 析 定 位 保护 限 值 大 于 告警 限 值 
但 定位 误差 小 于 定位 保护 限 值 的 概率 ， 不 考虑 漏 警 情况 : 
N(I(POSerraj| < HPL & HPL > HAL) 
Nan l 


上 面 各 式 中 ， 表示 满足 条 件 的 样本 个 数 ，Nn 表示 全 部 统计 样本 个 数 。 


(13) 


(14) 


HMIy = (15) 


Hunavai = (16) 


5 SISA 参数 完好 性 评 佑 试验 分 析 


本 文 利 用 BDS-3 系统 的 基本 完好 性 SISA 电文 和 IGMAS 全 球 站 观测 数据 ， 对 BDS-3 
系统 的 NPA 阶段 完好 性 服务 能 力 进行 分 析 ， 验 证 基本 完好 性 S154 电文 的 完好 性 性 能 。 根 
据 NPA 完好 性 服务 指标 要 求 ， 按 照 斯 坦 福 图 统计 方法 ， 评 估 BDS-3 系统 NPA 完好 性 服务 
可 用 性 及 完好 性 服务 风险 。 

利用 IGMAS 全 球 分 布 的 11 个 测 站 北斗 观测 数据 ， 系 统 完好 性 服务 性 能 评估 时 间 段 为 
2018 年 12 月 至 2020 年 10 月 ，IGMAS 测 站 分 布 如 图 3 所 示 。 

按照 BDS-3 运行 情况 ， 将 试验 验证 划分 为 基本 系统 阶段 和 完整 系统 两 个 阶段 ， 计 算 
IGMAS 站 的 水 平 保护 限 值 和 水 平定 位 误差 ， 分 析 BDS-3 系统 不 同 建设 阶段 NPA 完好 性 服 
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-135? 


评估 使 用 的 IGMAS 测 站 分 布 
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0° 


务 的 可 用 性 、 漏 警 率 和 严重 的 完好 性 风险 概率 。 在 试验 过 程 中 ， 没 有 出 现 “ 严 重 的 完好 性 风 


险 ” 事 件 ， 各 测 站 的 完好 性 服务 可 用 性 ( 即 包 络 能 


) 和 漏 警 率 统计 结果 见 表 4。 图 4 为 全 部 


测 站 NPA 完好 性 服务 可 用 性 的 综合 统计 结果 。 图 5 为 IGMAS 测 站 canb 2020 年 3 月 的 水 
平 误差 与 水 平 保护 限 值 时间 序 列 图 。 当 月 通过 SISA 参数 计算 出 的 五 PL 对 水 平 误差 的 包 络 


能 力 为 100%。 


表 4  IGMAS 各 测 站 的 完好 性 服务 可 用 性 ( 即 包 络 能 力 ) 和 漏 警 率 统计 


2019 年 7 月 一 12 月 


2020 年 7 月 一 11 月 


测 站 包 络 能 力 WERK | 测 站 包 络 能 漏 警 率 
/(%) /(%) /(%) /(%) 
Algr 98.77 0.000 | Byns 100.000 € 0.000 
Byns 99.81 0.027 | Brch 100.000 € 0.000 
Brch 99.58 0.000 | Canb 100.000 0.001 
Clgy 99.87 0.000 | Clgy 100.000 0.000 
Cnyr 99.97 0.000 | Cnyr 100.000 0.000 
Dwin 99.93 0.000 | Dwin 100.000 € 0.000 
Hmns 99.88 0.000 | Icuk 100.000 0.000 
Icuk 99.99 0.000 | Kndy 100.000 0.000 
Kndy 99.94 0.022 | Peth 100.000 0.000 
Peth 99.97 0.000 | Rdjn 100.000 0.000 
Taht 99.49 0.000 | Zhon 100.000 0.000 
平均 99.74 0.004 平均 100.000 0.000 
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J 日 存在 完好 性 
风险 的 概率 : 0% 


HPL/m 


0 非 精 密 进 近 服务 可 用 性 : 10096 7 f 


0 
误差 /m 
4 2020 年 7 月 至 11 月 NPA 完好 性 服务 可 用 性 统计 图 
40 
- Err,, Sta: canb 
| 包 络 能 力 : 100.00% |- 
: HPL HMITH-00000 . 
í 漏 警 率 ; 0.00% 
30r : z 误 警 率 : 0.0% 4 


水 平定 位 误差 与 水 平 保护 限 值 /m 


061 66 71 76 81 86 91 
2020 年 年 积 日 


图 5 2020 年 3 月 canb è HPL/HPE 时 间 序 列 图 


从 评估 结果 可 以 看 出 ， 随 着 星座 的 逐步 完善 ， 系 统 的 NPA 完好 性 服务 可 月 
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强 。 基 本 系统 阶段 BDS-3 系统 NPA 服务 可 用 性 为 99.75%， 全 星座 完整 系统 阶段 NPA 服务 
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可 用 


6 


广播 星 历 轨道 和 钟 差 信 息 的 预测 精度 。 由 于 系统 开通 服务 时 间 尚 短 ， 目 前 BDS-3 系统 对 用 


性 为 100.0096. 


结 论 


BDS-3 系统 采用 了 全 新 的 完好 性 设计 体制 ， 在 广播 星 历 中 ， 通 过 播发 SISA 参数 反映 


户 公开 的 ICD 文件 中 尚未 给 出 SISA 参数 的 用 户 算法 ， 影 响 了 系统 完好 性 服务 的 推广 应 用 。 


析 。 


BDS-3 系统 的 NPS 完好 性 服务 能 力 进行 了 初步 评估 。 实 验 结果 显示 ， 随 着 BDS-3 系统 星座 


本 文 设 计 了 SISAI 的 等 级 转换 方法 ， 并 研究 了 SISA 参数 的 用 户 应 用 算法 。 利 用 IGMAS 
全 球 观测 站 实测 数据 ， 采 用 斯 坦 福 图 分 析 方 法 ， 对 SISA 参数 的 完好 性 性 能 进行 了 验证 分 


根据 BDS-3 系统 的 完好 性 服务 目标 ， 按 照 国际 民航 NPA 阶段 完好 性 服务 指标 要 求 ， 对 


建设 的 完善 ， 系 统 的 完好 性 服务 能 力 增强 ， 在 全 部 星座 部 署 完 成 后 ，BDS-3 系统 对 NPA 阶 
段 的 完好 性 服务 能 力 为 100%。 本 文 研究 的 算法 和 试验 结果 可 为 深入 开展 BDS-3 系统 完好 性 


服务 


能 力 评估 提供 重要 参考 。 
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The User's Algorithm and Performance Evaluation on SISA 
of BDS-3 
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(1. Shanghai Astronomical Observation, Chinese Academy of Science, Shanghai 200030, China; 2. 
University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: The Beidou-3 global satellite navigation system was officially opened on the 
morning of July 31, 2020. With the gradual improvement of the Beidou system and the 
continuous improvement of its accuracy, users have higher and higher requirements for the 
integrity of the Beidou system. The BDS-3 system uses the four parameters (SISAI,., 
SISA», SISAl,4 and SISAI,;?) in the broadcast ephemeris to indicate the accuracy 
of the ephemeris and clock error parameters, but the four SIS AI parameters have not yet 
been given user algorithms in the interface control document (ICD). The integrity service 
of the BDS-3 system cannot be applied in practice. This paper studies the level conversion 
algorithm and user algorithm of the basic integrity parameter SISA of the Beidou-3 system, 
and uses the International GNSS Monitoring and Assessment System global network data 
and SIS AI parameters in the Beidou-3 satellite broadcast ephemeris in accordance with the 
integrity requirements of the International Civil Aviation Organization (ICAO) non-precision 
approach (NPA) and has conducted preliminary verification and analysis of the integrity and 
service performance of the BDS-3. The test results show that with the completion of the 
BDS-3 system and the SIS A parameter algorithm adopted in this paper, the BDS-3 system's 
NPA stage integrity service capability is 100%. 
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